# Metview Macro

t1000 = read("t1000.grb")

land_shade = mcoast(
    map_coastline_thickness         : 2,
    map_coastline_land_shade        : "on",
    map_coastline_land_shade_colour : "cream"
    )

polar_europe = geoview(
    map_projection      : "polar_stereographic",
    map_area_definition : "corners",
    area                : [30,-25,50,65],
    coastlines          : land_shade
    )

fixed_t = mcont(
    legend                         : "on",
    contour                        : 'off',
    contour_level_selection_type   : "level_list",
    contour_level_list             : [-50,-20,-5,0,5,10,15,20,25,50],
    contour_shade                  : "on",
    contour_shade_method           : "area_fill",
    contour_shade_max_level_colour : "red",
    contour_shade_min_level_colour : "blue",
    contour_shade_colour_direction : "clockwise"
    )

plot(polar_europe, t1000, fixed_t)